Signature of electronic excitations in the Raman spectrum of graphene. 
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Inelastic light scattering from Dirac-type electrons in graphene is shown to be dominated by the 
generation of the inter-band electronic modes which are odd in terms of time-inversion symmetry 
and belong to the irreducible representation A2 of the point group C6« of the honeycomb crystal. At 
high magnetic fields, these electron-hole excitations appear as peculiar n~ — > n + inter-Landau-level 
modes with energies Lo n — 2\/2n hv/Xs and characteristically crossed polarisation of in/out photons. 

PACS numbers: 73.63.Bd, 71.70.Di, 73.43.Cd, 81.05.Uw 



O 

(N 
C 

in 

(N 



1 

c 

o 
o 



> 
in 

(N 

in 

\6 
o 

OS 

o 



X 



Inelastic (Raman) scattering of light is a powerful tool 
to study excitations in solids Recently, Raman spec- 
troscopy has been used to study phonons in graphene Q , 
where it has become the method of choice for determin- 
ing the number of atomic layers in graphitic flakes 043 ■ 
In particular, single- and multiple-phonon-emission lines 
in the Raman spectrum of graphene and the influence 
of the electron-phonon coupling on the phonon spectrum 
have been investigated in great detail [9HL6I]. However, no 
experimental observation or comprehensive theoretical 
analysis has been reported on the Raman spectroscopy 
of electronic excitations in graphene, despite extensive 
studies of optical absorption in this material fl7l - l23l |. 



In this Letter, we present a theory of inelastic light 
scattering in the visible range of photon energies accom- 
panied by electronic excitations in graphene. We classify 
the relevant modes according to their symmetry and pre- 
dict peculiar selection rules for the Raman-active excita- 
tions of electrons between Landau levels in graphene at 
quantizing magnetic fields. Graphene is a gapless semi- 
conductor [U, [25[ , with an almost linear Dirac-type spec- 
trum, £ = avp in the conduction (a = +) and valence 
(a = — ) band, which touch each other in the corners of 
the hexagonal Brillouin zone, usually called valleys. The 
band structure of graphene is prescribed by the hexago- 
nal symmetry Cq v of its honeycomb lattice, and it is nat- 
ural to relate Raman-active modes to the irreducible rep- 
resentations of the point group C6„. We argue that the 
dominant electronic modes generated by inelastic scat- 
tering of photons with energy less than the bandwidth 
of graphene are superpositions of the interband electron- 
hole pairs which have symmetry of the representation 
A2 of the group Cq v and are odd with respect to the 
inversion of time. Their excitation process consists of 
two steps: the absorption (emission) of a photon with 
energy Q (Q = £1 — lS) transfering an electron from an 
occupied state in the valence band into a virtual state 
in the conduction band, followed by emission (absorp- 
tion) of the second photon with energy £1 (f2). Its am- 
plitude is determined by the sum of partial amplitudes 
distinguished by the order of absorption and emission of 
photons, and by which carrier in the intermediate state 
(an electron above the Fermi level or hole below it) un- 
dergoes the second optical transition. The dominance of 
such process over the process involving the contact in- 




FIG. 1: Spectral density g{oS) of light inelastically scatterred 
from electronic excitations in graphene at quantising magnetic 
fields (solid line) and at B — (dashed line). Here ui <C O. 
is the Raman shift. Sketch illustrates intermediate and final 
states of the dominant Raman process. 



teraction [26|, [27( of an electron with two photons is a 
peculiarity of the Dirac-type electrons in graphene. Fill- 
ing the conduction band or depleting the valence band, 
up to the Fermi level a[i, forbids the excitation of inter- 
band electron-hole pairs with energies to < 2/j, leading to 
the linear Raman spectrum with a 2 /i threshold [dashed 
line in Fig. [T] . The quantization of the electronic spec- 
trum into Landau levels [3(| (LL) e[n a ] = a\/2nhv/\B 
in a strong magnetic field (As = ^Jhc/eB is the mag- 
netic length, n — 0, 1, 2, ...) makes Raman spectrum dis- 
crete at low energies uj n = 2^/2nhv / Xb, with peculiar for 
the Dirac-type electrons selection rules, rT — > n + of the 
dominant Raman-active transitions [solid line in Fig. [T]. 

The following theory is based upon the tight-binding 
model of electron states in graphene expanded into the 
Dirac-type Hamiltonian |3lj . 

v 2 

n = v^-P A^EP^SP)^. (1) 
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The latter describes electrons in the conduction and 
valence bands around the Brillouin zone (BZ) corners 
K and K'. We use notations [321 such that S = 



(P'kk 

a z = n 
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where a 



x/y/z 



'AB 



and 



and lfj( v J, z are Pauli 



KK' w ± AB, wiicic u AB 

matrices acting on A-B (sublattice) and K-K 1 (valley) 
indices of the four-component wave function {ipK, ipK'}, 
where ipK — Vpa,<Pb] and ipK' — Vpb,^a]- While 4- 
spinors {ipK^K 1 } realise 4D irreducible representation 
of the full symmetry group of the crystal, the valley- 
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TABLE I: Cq v irreducible representations by the valley- 
diagonal operators E 8 and A Z E\ 



diagonal operators Yi 1 and A Z E Z can be combined into 
irreducible representations d, HH of the group Ce^, Ta- 
ble HI and A z is used to describe valley- asymmetry of 
Dirac electrons. The first term in T-L determines the lin- 
ear spectrum avp with v ss 10 8 cm/s and p being the 
in-plane momentum counted from the BZ corner. The 
second term takes into account weak trigonal warping 
[hopping parameter 70 ~ 3eV determines the bandwidth, 
~ 670], which has an inverted shape in the opposite 
corners of the BZ [31]. The vector potential of light 
A = Ei. q . ?I (le i(qr_nt)/R &q, g „i + h.c.) is included in 
P = p -A, where & q ,g z ,i annihilates a photon char- 
acterised by the polarisation [1 for incident and 1 for 
scattered light], in-plane momentum q, energy f2, and 
q z = ^Q 2 /c 2 -q 2 . 

The amplitude R = Rd + R w + TV of the Raman 
process with the excitation of an electron-hole (e-h) pair 
in the final state corresponds to the Feynman diagrams 
shown in Fig. [5J Here, we call an 'electron' an excited 
quasiparticle above the Fermi level aji, and a 'hole' an 
empty state at e < afi. The building blocks of the dia- 
grams include Green's functions for the electrons and the 
electron-photon interaction vertices: 
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In the amplitude i?, the term represents the contri- 
bution of the first two diagrams in Fig. [2] They describe 
a photon- assisted transition of an electron with momen- 
tum p from under the Fermi level into a strongly off- 
resonant virtual intermediate state (note that w|p + k| ss 
vp <C f2), followed by another transition (of either 

electron or a hole) which returns the system onto the en- 
ergy shell. The two diagrams in Rd differ by the order 
of absorption/emission of the photons with J7,f2 3> vp, 
and, therefore, by the sign of the energy denominator 
in G R I A . A partial cancellation between them deter- 
mines the effective 2-photon coupling to the electrons 
characterised by the matrix form in the representation 
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FIG. 2: Feynman diagrams describing Raman scattering with 
the excitations of electron-hole pairs in the final state. 



A2, Table [H As a result, such process excites a 'valley- 
symmetric' electronic mode corresponding to the repre- 
sentation A2 of C6„ and odd in terms of time-inversion 
symmetry. The term R w in Eq. describes the contact 
interaction between an electron and two photons charac- 
terised by d 2 R/dpidpj. Although for free non-relativistic 
electrons contact interaction is important 26] , for Dirac- 
type electrons it is absent. It reappears only after devi- 
ations from the Dirac spectrum are taken into account, 
i.e., the 'valley-antisymmetric' warping term in Eq. <jTJ) , 
and R w generates excitations with the symmetry of the 
representation E2 of Cq v . For scattering of photons with 
Q < 70, Rw < Rd [27]. Finally, TV stands for the con- 
tribution of the diagrams containing a 'triangular' loop 
T and the RPA-screened electron-electron interaction 
V. It accounts for the generation of a virtual e-h pair 
which recombincs creating a real c-h excitation through 
the electron-electron interaction, and its effect is negligi- 
bly small [H|. 

The probability for a photon to undergo inelastic scat- 
tering from the state (q, q z ) with energy f2 into a state 
(q = q — k, q z ) with energy Q = — cj, by exciting an 
e-h pair in graphene with Fermi energy afx at low tem- 
perature T < lj, is 



d 2 p 



/p"«(l - /(p+k)" - £ (p+k)° +W) 



Airh 3 ' 

x J2 tr + V^n p )R + (l + aSn p+k )} . (3) 

Here, a = ± distinguishes between n- and p-doping 
of graphene, 77 = — /+ stands for the excitation of the 
inter/intra-band electron- hole pairs, and spin-degeneracy 
is taken into account. The probability w describes the 
angle-resolved Raman spectrum, as opposed to the angle- 
integrated spectral density, 
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In undoped graphene the inter-band e-h pairs are the 
only allowed electronic excitations. The probability, 
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of their excitation by photons with f2 < 70 is dominated 
by the contribution, Rjj of the first two diagrams in 
Fig. [2] This determines typically crossed linear polar- 
isation of in/out photons described by the polarisation 
factor S s , which is equivalent to saying that they have the 
same circular polarisation, in contrast to a weak contribu- 
tion of the process enabled by the warping term (second 
term in Eq. (|5]l). with the opposite circular polarisation 
of in and out photons described by the factor S . 

In doped graphene, with fi 3> Clv/c, inter-band elec- 
tronic excitations with uj < 2\x are blocked, so that 
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After integrating over all directions of the propagation 
of scattered photons, we find the spectral density of the 
angle-integrated Raman signal, 
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Here F(\x\ < 1) = |(1 + x) and F(\x\ > 1) = 9{x), 
step function. In undoped graphene (fi = 0), spectral 
density g(ui) corresponds to the yield 1$ = g(oj)doj ~ 

(tt c I fy) such that ~ ~ ltr 10 . 

In doped graphene one may also expect to see some 
manifestation of the intra-band e-h excitations in the 
vicinity of the Fermi level, with a small energy transfer 
u) < Qv/c. Their analysis requires taking into account 
all diagrams in Fig. [21 in particular, due to an additional 
asymmetry between the conduction and valence bands 
caused by the difference of their filling which increases 
the value of the triangular loop 34| , 
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Then, we find that, for ui < (w/c)0 <C < 70, 
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The yield of this low-energy feature is SI = J Sg(uj)duj ~ 
10" 15 for fl ~ leV |36|. 

Electronic spectrum of graphene in a strong magnetic 
field can be described as a sequence n a of Landau levels 
(LLs). sjn a ] = ae n with e n = \/2nhv/XB, corresponding 
to [U |3[J] the states \n a ) = ^($„, ia$„_i) for n > 1 

and |0) = ($o,0) (where Xb = \JhcjeB and $„ are 
the normalised LL wave functions in the Landau gauge). 



Then, electron's Green functions and interaction vertices 
leading to optically active inter-LL excitations in mono- 
layer graphene summarised in Table [TT] take the form 
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Here e± = (e^ ± ie y ) is used to stress that a circularly 
polarised photon carries angular momentum m = ±1. 

The excitation of the e-h pairs by Raman scattering 
in graphene at strong magnetic fields characterised by 
the first two Feynmann diagrams in Fig. [5] produces the 
electronic transition n~ — > n + between LLs, with angular 
momentum transfer Am = and excitation energy lu = 
2e n [Fig. [I], and transitions (n — 1)~ — > [n + 1) + and 
(n+l)~ — > (n— 1) + , with Am = ±2 and u; = £„_i+£„ + i. 
The amplitudes of these two processes, 
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are such that R n -^ n + ^ -R( n =pi)--i.(n±i)+ for cj O, 
due to a partial cancellation of the two diagrams consti- 
tuting Rjj. Notice that these inter-LL modes n~ —> n + 
have the symmetry of the representation A2 in Table U 
and the same circular polarisation of in and out photons 
involved in its excitation. Finally, the contact term R w 
in Fig. fallows for a weak transition n" — > (n± 1) + , with 
the amplitude R w <C R n -^ n + Superficially, such a 
transition, with Am = ±1 resembles the inter-LL transi- 
tion involved in the far-infrared (FIR) absorption |17ll2l| . 
However, the FIR-active excitation is 'valley-symmetric' 
[2~D ] and corresponds to the representation Ei, whereas 
the Raman- active rf~ — > (n ± 1) + mode corresponds to 
E2, allowing the latter to couple to the T-point optical 
phonon and, thus, leading to the magneto-phonon res- 
onance feature in the Raman spectrum Also, R w 
originates from the trigonal warping term in % which vi- 
olates the rotational symmetry of the Dirac Hamiltonian 
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weak in Raman, 
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A2 n~ -> n+ dominant in Raman a ± — > it 



TABLE II: Raman-active inter-LL excitations in graphene. 

by transferring angular momentum ±3 from electrons to 
the lattice, so that initial and final state photons in it 
have opposite circular polarisations. 

The dominant inter-LL transitions n" — > n + deter- 
mines the spectral density of light scattered from elec- 
tronic excitations in graphene at high magnetic fields: 

ffn-->n+M ~ 5 S (X 6 Iq 8 ^) ^ln{u-U n ). (8) 
^ C ^ ' n>l 

Here j n (x) = Tr^ 1 T n /[x 2 + r 2 ], and T n is inelastic LL 



broadening which increases with the LL number, u> n — 
2e n = 2V2nfrv/\ B , and the factor E s = |lxl*| 2 in Eq. © 
indicates that in and out photons have the same circular 
polarisation. 

The "nT — > n + inter-LL transitions are specific for 
Dirac-type electrons in graphene and represent the most 
pronounced signature of its electronic excitations in the 
Raman spectrum. The quantum efficiency of the low- 
est, u\ — 2^/2Hv/\b peak in the spectrum in Fig. [1] is 

J i ~ ^-^r) P er incoming photon. For B = 20T, 
we estimate I\ ~ 10~ 12 for photons with energies in the 
visible range, which is feasible to detect in the inelastic 
light scattering experiments. 
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